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Abstract 

The  broad  goal  of  our  research  program  was  to  study  excited  state  chem¬ 
istry,  and  electronic  energy  transfer  and  chemical  energy  redistribution 
processes,  using  the  methods  of  ab  initio  electronic  structure  theory.  In 
the  course  of  this  research  we  have  contributed  to  the  rapid  growth  in  muiti- 
configuration_^if-^nsistent_^ield  (MCSCF)  methodology,  initially,  with  the 
development  of  an  MCSCF  procedure  based  on  the  generalized  Brillouins  theorem 
and  the  approximate  natural  orbital  method  (GBT-N'O  MCSCF)  and  more  recently 
with  the  extension  of  Lengsfield's  density  matrix  driven  quadratically  con¬ 
vergent  MCSCF  procedure  to  treat  state  averaged  wavefunctions  (SA  MCSCF) . 

The  later  algorithm  has  significantly  expanded  the  scope  of  problems  in 
excited  state  chemistry  amenable  to  treatment  at  the  MCSCF/CI  level. 

As  a  prototypical  system  we  chose  reactions  of  the  form 

Me(Y)  +  N'  0(4“^)  ^  MeO(Z)  +  N'.,(^Z'^)  I 

2  2  g 
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where  Me(\')  =  Mg(  3,  P)  and  Ca(  S,  ?) .  Although  both  the  macnesium  ana 
calcium  reactions  are  highly  exoerglc  the  reactivity  and  the  distributicn 
of  product  electronic  states  differs  narkedl;.-  for  these  two  homologues. 

Over  the  period  of  this  grant  we  have  studied  the  properties  of  several 
subspecies  participating  in  this  reaction  including  raaiatlve  decav  or 
metastable  atomic  species,  and  the  electronic  structure,  bonding  and  optica, 
properties  of  the  alkaline  earth  oxides.  L’ltimately  we  were  able  to  characterize 
several  regions  of  the  lowest  singlet  and  triplet  poten  Liul  cuor-v  surrucii.s 
of  reaction  I  at  che  SC^/ two  configuration  SCFCICSCI- ) /Cl  level. 

.1.  Radiative  Properties  of  Ca(^D,  '^?) 

Using  GBT-N’O  MCSCF/CI  metlrods  tiie  Cs,  1^’^P  and  1^'^D  states  of  calcium 
atom  were  characterized.  Using  an  ort.hogonal  orbital  transition  moment 
method  and  atomic  spectral  data,  line  strengths  for  ~  Ip  -  Ig,  'P  *  *3 
and  other  transitions  were  estimated.  It  was  shown  that  che  principal  decav 
channel  of  Cal^DT)  is  Che  spin  forbidden  dipole  allowed  path  through  che 
state  to  che  state. 


B.  Electronic  Structure  of  .Alkaline  Earth  Oxides 

MCSCF/CI  wavefunctions  were  determined  for  BeO,  MgO  and  CaO  with  partic¬ 
ular  attention  focused  on  che  1,2  states.  It  was  shotvn  that  a  viable 
starting  point  for  the  description  of  these  states  is  necessarily  multicon- 
figurational  in  character.  In  each  case  contributions  from  structures  .Me~-0“2 
and  Me'*'0“  must  be  considered  simultaneously,  and  in  the  case  of  the  Me‘*'0~ 
structures  alignment  of  the  oxygen  p-hole  parallel  and  perpendicular  to  the 
axis  must  be  included.  Subsequently  dipole  allowed  radiative  transitions 
were  studied  for  the  2  and  i  2ri  states  in  MgO  and  CaO.  Here  the  SA-MCSCF 
procedure  oroved  to  be  essential. 
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C.  The  Reaction  Mg  +  N^O  =  MgO  + 


The  lowest  surfaces  of  this  reaction  were  characterized  usinc 

SCF/TCSCF/CI  methods.  These  calculations,  which  were  made  economically 
feasible  with  the  purchase  (with  funds  from  AFOSR  and  N'SF)  of  a  dedicated 
minicomputer  system,  have  enabled  us  to  characterize  the  mechanism  for  the 
production  of  MgO(X  and  a  3n)  and  also  explain  the  disparities  between 
the  magnesium  and  calcium  reactions.  The  importance  of  curve  crossings  with 
open  shell  structures  corresponding  to  excited  atomic  states  in  the  grounc 
state  reaction  was  demonstrated.  A  clear  phenomenological  description  of 
this  reaction  has  emerged.  The  ^’'^A*  surfaces  can  be  qualitatively  parti¬ 
tioned  into  three  contiguous,  nondisjoint  regions  characterized  by  distinct 
values  of  an  approximate  reaction  coordinate,  a  reactant  rezion  in  whicii 
s  =  -  R(Mg-O) ,  a  bending  reaction  in  which  ;  =  and  a  product 


region  in  which 


R(N-O).  Followin; 


are  driven  by  a  Rener-Teller  instability  in  M^0~(^r.). 
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Project  Summary 

Objectives  and  Accomplishments 

The  understanding  of  the  properties  and  reactivity  of  transient  elec¬ 
tronically  excited  chemical  species  has  important  implications  for  the  devel¬ 
opment  of  chemical  laser  systems,  and  energy  storage  and  detection  devices. 
Our  approach  to  the  study  of  excited  state  chemistry  is  theoretical  using 
the  methods  of  ab  initio  electronic  structure  theory,  in  particular  stated 
averaged  multiconfiguration  self-consistent  field  (SA-MCSCF)  and  configura¬ 
tion  interaction  (Cl)  methods  to  characterize  the  relevant  wavef unctions  with 
the  Born  Oppenheimer  approximation. 

(A)  Computational  Methods 

In  the  treatment  of  the  reactive  systems  involving  the  interaction  of 
several  diabatic  surfaces  or  excited  electronic  states  the  MCSCF  approxima¬ 
tion  provides  a  valuable  starting  point,  Unfortuneately  the  range  of  prob¬ 
lems  in  excited  state  chemistry  amenable  to  treatment  at  the  MCSCF  level  is 
limited  by  variational  collapse  which  precludes  direct  optimization  of  ex¬ 
cited  states  of  given  electronic  symmetry.  Our  studies  of  the  SA-MCSCF  pro¬ 
cedure  (bibliography  references  (BR)  7,8)  have  shown  this  approach  to  provide 
a  generally  useful  method  for  circumventing  this  difficulty  thereby  extending 
the  scope  of  MCSCF  methodology. 

To  this  point  we  have  concentrated  on  the  use  of  this  procedure  in  the 
evaluation  of  interstate  properties,  total  energies,  spectroscopic  constants 
and  dipole  moments  and  nonrelativistic  interstate  matrix  elements,  in  par¬ 
ticular  transition  dipole  moments.  The  later  have,  in  turn,  been  used  to 
study  optical  properties  of  the  alkaline  eath  oxides  (BR  7-9) .  However  a 
potentially  significant  application  of  this  approach  exists  in  the  area  of 
spin  forbidden  chemistry  in  which  SA-MCSCF/CI  methods  could  be  applied  to 
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the  determination  of  interstate  matrix  elements  of  the  spin  orbit  operator. 

An  application  of  this  approach  to  the  reactions 

Mg(^P)  +  Ca(^S)  =  Ca(^D)  +  Mg(^S) 

=  Ca(^D)  +  Mg(^S) 

is  discussed  in  a  proposal  recently  submitted  to  AFOSR. 

(B)  Dedicated  Minicomputer  System 

Ab  initio  electronic  structure  methods  require  large  amounts  of  computer 
time.  In  the  absence  of  a  (subsidized)  readily  accessible  mainframe  computer 
system  reserved  for  various  forms  of  scientific  computation  a  dedicated  mini¬ 
computer  system  (i.e.  a  computer  dedicated  to  a  single  class  or  limited  number 
of  related  classes  of  problems)  represents  an  economically  viable  approach 
to  obtaining  the  requisite  computational  facilities.  Here  the  economies 
result  principally  from  reduced  overhead  associated  a  user  run  system  and 
the  increased  operating  efficiency  resulting  from  improved  scheduling. 

During  the  course  of  this  grant  we  have  been  fortunate  to  acquire  (prin¬ 
cipally  with  funds  from  AFOSR  and  NSF)  a  Perkin  Elmer  3200  series  computer 
system,  presently  configured  with  a  PE  3230,  1.5  megabyte  (MB)  central  pro¬ 
cessing  unit  (CPU),  2-80  MB  disc  drives,  1600  bit/inch  (BPI)  tape  drive  and 
a  line  printer.  Although  as  discussed  in  a  University  Research  Equipment 
Proposal,  additional  hardware  would  significantly  increase  the 
scope  of  problems  accessible  to  this  system,  this  system  despite  its  recent 
acquisition,  has  been  indispensible  in  achieving  the  goals  of  our  research 
program.  In  particular  the  characterization  of  the  reaction  (BR  10) 

M8(^S,^P)  +  N„0(^Z‘^)  =  Mg0(.X^i:'^,a^H)  +  N_(^:^) 
discussed  below,  would  have  been  (economically)  impossible  without  tlie  PE 
system. 
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(C)  Excited  State  Chemistry 

The  chemistry  studied  during  the  course  of  this  grant  fails  naturallv  into 
three  categories  which  are  summarized  below: 


a.  Radiative  Decay  of  Metastable  States  of  Ca  Spin  Forbidden  and  Spin 
Allowed  Mechanisms 

Several  experimental  investigations  have  uncovered  apparent  anamolies 

in  radiative  decay  processes  involving  excited  states  of  Ca,  Sr  and  Ba.  ^s 

a  first  attempt  to  understand  these  processes  we  used  our  Generalized 

Brilburns  theorem  (GBT)  MCSCF/CI  methodology,  together  with  accurate  atomic 

spectra  data  to  determine  branching  ratios  for  radiative  transitions  out  of 
3  3  1  1 

the  lowest  P,  D,  D  and  P  states  of  Ca  (BR  3).  Both  spin-forbidden  dipole- 
allowed  (SFDA)  and  spin-allowed  dipole-forbidden  (SADF)  mechanisms  were  con¬ 
sidered  in  addition  to  the  usual  dipole-  and  spin-allowed  channel.  Relativ¬ 
istic  effects  were  included  in  an  approximate  seroi-erapirical  manner.  The 
results  of  this  investigation  were  quite  encouraging.  Tlie  decay  of  metastable 
Ca(^D)  was  shown  to  proceed  largely  by  a  SFDA  mechanism  to  the  Ca(^?)  state 
rather  then  as  might  have  been  expected  by  direct  quadrupolar  decay  (a  S.VDF 
mechanism)  to  the  ground  state.  In  addition  these  results  suggest  a  pos¬ 
sible  explanation  for  the  weak  lasting  action  of  the  (^P-^D)  compared  to  the 
3  3 

strong  (  D-  P)-transition  in  Sr  vapor. 

b.  Electronic  Structure  and  Radiative  Properties  of  Alkaline  Earth  Oxides 
A  series  of  studies  of  the  low-lying  states  of  BeO,  MgO  and  CaO  was 
performed  with  particular  emphasis  on  the  1,2,3  states  for  which  treat¬ 
ment  at  the  (SA)-MCSCF/CI  level  appears  to  be  essential.  It  luis  emerged  that 


3 


r  - 


the  electronic  structure  of  these  homologous  states  is  qualitatively  distinct. 

The  differences  were  found  to  result  from  a  combination  of  three  factors,  the 

+  +2  -2 
stability  of  the  states  of  Me  ,  the  relative  contributions  from  Me  0  and 

Me'^'^O  ^  structures  and  the  orientation  of  the  oxygen  p-hole,  or  p..,  in  the 

+1  -1 

Me  0  structure.  Spectroscopic  constants  for  the  as  yet  unobserved  b 
state  of  MgO  were  determined.  In  particular 

(i)  BeO  (BR  1,4) 

Using  a  Gaussian  basis  set  of  slightly  better  than  double  zeta  plus 
polarization  quality,  MCSCF  and  Cl  calculations  were  performed  for  the 
ground  and  several  low-lying  singlet  states  of  BeO.  The  computed  sepa¬ 
rations  for  A^n  and  B^I!^  states  are  in  good  agreement  with  experi¬ 

ment.  The  and  states,  the  likely  upper  states  in  the  Bengtsson- 
Knave  and  Harvey-Bell  systems,  are  predicted  to  lie  '''40,000  cm”^  above 
the  X  state,  while  the  2^11  and  3^1"*"  are  predicted  to  lie  56,000  and 
65,000  cm  ^  above  the  X  state,  respectively. 

(ii)  MgO  (BR  2, 5, 7, 9) 

Using  a  double  zeta  plus  polarization  basis  set  of  Slater  orbitals,  full 
valence  MCSCF  (FVMCSCF)  calculations  were  performed  for  the  low-lying  states 
of  MgO.  For  each  state  FVMCSCF  calculations  were  used  to  identify  the  important 
configurations  which  were  then  used  in  an  MCSCF  calculation  and  subse¬ 
quently  as  references  in  a  single  and  double  excitation  Cl  calculation. 

This  approach  was  found  to  treat  all  states  equivalently,  with  the  maximum 

error  in  the  computed  T^'s  and  of  800  cm  ^  and  '^0.03.^,  resoectively . 

3  + 

The  b  Z  state  which  has  yet  to  be  characterized  experimentally  vas  pre¬ 
dicted  to  have  a  T^  of  '8300  cm  ^  and  a  bond  length  of  1.79  A.  A  spec- 


Lroscopii:  inalvsis  iM’  the  potential  curves  indicates  that  their  shapes 


are  in  quite  reasonable  agreement  wi^’  the  range  of  experimental  results. 
Dipole  and  transition  dipole  moments  were  determined  for  the  X,B  and 
states  at  the  MCSCF/CI  level.  Wavefunction  determined  at  the  SA- 
MCSCF/CI  level  were  used  to  discuss  radiative  decay  of  the  and 

A^n  states,  i.e.  the  B^I  “^AH,  BZ  '’•X_  and  A  .1  X  radiative 
transitions. 

(iii)  CaO  (BR  8) 

State  averaged  multiconfiguration  self-consistent  field  and  configura¬ 
tion  interaction  methods  were  used  to  discuss  the  1,2^1^  states  of  CaO. 

Spectroscopic  constants  R  ,  T^,  and  were  determined .  In  addition, 

1  +  1 

radiative  transitions  between  the  1,2  Z  states,  and  the  1  x  states 

were  discussed.  The  spectroscopic  constants  are  in  reasonable  agreement 

with  experiment  although  the  predicted  bond  lengths  are  uniformly  long. 

1  + 

The  calculated  lifetime  for  the  v  =  6  level  of  the  A  I  state  of  130  ■ 

-9 

10  s  is  in  good  agreement  with  a  preliminary  experimental  value.  The 
lifetime  of  low-lying  vibrational  states  in  the  A' manifold  are  over 
an  order  of  magnitude  longer  than  their  A  state  counterparts.  The  = 

0  level  of  the  A'  state  is  predicted  to  have  a  lifetime  of  25  •  10  s 


c.  The  Reaction  Me  +■  =  MeO  + 

^  2  g 

The  lowest  singlet  and  triplet  potential  energy  surfaces  for  this  reac¬ 
tion  have  been  investigated  for  Me  =  Mg(BR  10).  Several  regions  of  the  C 

s 

11  3  3 

symmetry  surfaces  corresponding  to  Mg(  S)[l  A']  and  Mg(^P)[l  A']  were  con- 

3  3 

sidered.  On  the  1  surface,  which  correlates  diabatically  with  Mg0(a  7.), 

the  reaction  occurs  via  charge  transfer  from  magnesium  to  colinear  N'  0,  i.e. 
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Mg  +  N^O  becomes  .N',,0  -  Mg  .  Kquivaiently  the  highest  singly  occupied 

orbital  changes  from  magnesium  dp_.  to  M,,U  J:  in  character.  in  the  -actant 

channel  r-type  approach  is  doimhill  and  favored  over  the  alternative  uphill 
J-type  approach.  Charge  transfer  leads  to  distortion  of  from  its  equil¬ 

ibrium  geometry  to  a  bent  structure.  Entrance  into  the  product  channel 
then  corresponds  largely  to  motion  along  the  X-0  coordinate  and  involves  an 

orbital  reorganization  in  the  N^O  moiety.  Equivalently  this  reorganization 

3  3 

represents  an  avoided  crossing  of  A*  surfaces  correlating  with  MgO(a  + 


and  .MgOCb^l"^)  +  In  this  region,  which  is  stable  with  re- 


spect  to  the  reactant  asymptote  R(N-N)  is  larger  than  in  N'  (  Z  ). 

-  S 

On  the  1  a'  surface,  there  is  a  barrier  to  charge  transfer.  .Approach  of 
Mg(^S)  to  N'20  is  uphill.  The  barrier  to  charge  transfer  is  associated  with 
the  absence  of  accessible  open  shell  configurations  which  correlate  with  e.x- 
cited  singlet  states  in  the  isolated  atom.  For  all  nuclear  configurations 


sampled  the  l\\'  surface  is  lower  than  tiie  surface.  However  followin. 

1  3 

charge  transfer  the  separation  of  the  1  .\'  and  1  A'  surface  decreases  marh- 
edlv. 


3,. 


In  the  Ca+-  N^O  system,  Ca(  S)  is  reactive.  For  calcium  the  e.xcited 
singlet  states  have  significantly  lower  excitation  energies  than  their  counter 
parts  in  magnesium.  This  suggests  that  the  mechanism  for  the  reaction  of 
Ca(^S)  involves  a  preliminary  crossing  with  surfaces  correlating  with  Ca(^D) 

1  3 

or  Ca(  P)  followed  by  the  charge  transfer  mechanism  discussed  above  for  Mg(  •’) 


Single  configuration  (non  variational)  SCF  calculations  suggest  that  the 
electronic  path  leading  to  the  Mgi)  arc  band  states  is  uphill,  a  feature 


which  may  be  attributable  to  lar^e  energy  difference  between  Mg  (  S)  and 

+  2  1 

Mg  (  f) .  Tlie  Ca  (  r)  +  N^O  reaction  is  known  to  produce  chemiluminescence 
attributable  to  arc  band  emission.  In  calcium,  the  Ca^(^S)  -  Ca^(^P)  separa¬ 
tion  is  signif icantlv  (-10000  cm  S  smaller  suggesting  an  additional  direction 
for  future  theoretical  investigation. 
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